Gap soliton formation by nonlinear supratransmission in Bragg media 
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A Bragg medium in the nonlinear Kerr regime, submitted to incident cw-radiation at a frequency 
in a band gap, switches from total reflection to transmission when the incident energy overcomes 
some threshold. We demonstrate that this is a result of nonlinear supratransmission, which allows 
to prove that : i) the threshold incident amplitude is simply expressed in terms of the deviation from 
the Bragg resonance, ii) the process is not the result of a shift of the gap in the nonlinear dispersion 
relation, iii) the transmission does occur by means of gap soliton trains, as experimentally observed 
[D. Taverner et at, Opt Lett 23 (1998) 328], iv) the required energy tends to zero close to the band 
edge. 



Introduction. Light propagation in dielectric media 
with periodically varying index acquires intriguing prop- 
erties when Kerr nonlinearities come into play. From a 
linear theory H we know that such a medium is a pho- 
tonic band gap structure, or Bragg medium, which is to- 
tally reflective when the frequency of the incoming wave 
falls in one of the gaps. 

However, for sufficiently energetic irradiation, the Kerr 
effect becomes sensible and the Bragg medium may be- 
come transparent in the gap, switching from total reflec- 
tion to high transmissivity, which has been predicted in 
1979 |2| by using the stationary coupled mode approach 
of , and which has been experimentally realized for the 
first time in 1992 [i]. 

In 1987 the existence of effective wave propagation in 
a gap has been associated to the existence of gap soli- 
tons derived from slowly varying envelope limits of the 
Maxwell equation in a periodic structure 0, |^. For 
unidirectionnal propagation, fast Kerr nonlinearity and 
slowly varying evelope approximation, one of these limit 
models is the coupled mode system (Tj] which describes 
the evolution of the envelopes of the contrapropagating 
electric fields. A rigorous derivation of from the anhar- 
monic Maxwell-Lorentz equations by the method of mul- 
tiple scales is given in Q. The explicit soliton-like solu- 
tion of the coupled mode equation has been discovered 
in 1989 0, and at the same time generalized to a two- 
parameter family in . This two-parameter gap soliton 
solution, that we use here, has motivated many experi- 
mental searches of the Bragg soliton. 

After the pioneering experiments of in steady-state 
regime, the first experimental observations of nonlinear 
pulse propagation in a fiber Bragg grating have been per- 
formed in 1996 'll] under laser pulse irradiation at a 
frequency near (but not inside) the photonic band gap. 
In that case the nonlinearity acts as a source of modula- 
tionnal instability generating trains of pulses T^] . Partic- 
ularly interesting is the observation of pulse propagation 
at a velocity less than the light velocity. Later in 1998 
|13| a convincing experiment with a fiber Bragg grating 



demonstrated the repeated formation of gap soliton un- 
der quasi-constant wave irradiation inside a band gap, 
this is the problem we are interested in. 

There has been a great deal of discussions concerning 
the exact role of the gap soliton in the switching to high 
transmissivity, see e.g. 14;|. A rather natural intuition 
is that the nonlinearity shifts the gap and hence allows 
the medium to become transparent and to create a gap 
soliton. We shall discover that this intuition is wrong as 
the process at the origin of high transmissivity in Bragg 
media is the nonlinear supratransmission recently discov- 
ered in the sine-Gordon chain p^ll6j and experimentally 
realized (together with application to Josephson junc- 
tions arrays) in 0. As a consequence, the switching to 
high transmissivity is indeed accomplished by means of 
gap soliton generation as a result of a fundamental insta- 
bility [13 which are effectively the objects experimentally 
detected in the output in p!^ . 

This result allows us to predict analytically by formula 
(ISJ the threshold of incident energy of a cw-radiation 
above which high transmissivity is reached, as an ex- 
plicit simple function of the departure (denoted by the 
dimensionless angular frequency Q) from the Bragg fre- 
quency (gap center). In particular we obain that the pro- 
cess requires less energy for frequencies close to the band 
edge f2 = 1, a result which was previously attributed to 
pulse shape properties. Interestingly enough, in the limit 
51 — !■ 1 the required energy vanishes, opening the way to 
experimental realizations. 

The model and its solution. Our starting point is the 
basic coupled mode system 0, Q governing the forward 
Ef and backward Eb slowly varying envelopes of the elec- 
tric field 
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E{Z, T) is the transverse polarized component propagat- 
ing in the direction Z at frequency close to the Bragg 
frequency ojq. Following j^] we write the coupled mode 
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system in reduced units as 

4+|i+/+(y+i/p).^o, 

.||-fl + . + (> + |eR/^» 



(2) 



The reduced units are z = kZ and t = kcT where c 
is the hght velocity in the medium and n the couphng 
constant. The reduced field variables are e = Efyj2T/K 
and / = Ei,^j2T / HI where F is the nonlinear factor. The 



(linear) dispersion relation of (|2Jl, namely 



1 + 



possess the gap [—1,+!]. 

The solitary wave solution of |(5J) given in has a 
simple stationnary expression 

e = y273A e-^^"*-^) cosh^^ [A(z - ^o) - ^g] , (3) 

i 



f = -V273Ae-^("*-'^) cosh-i[A(z - zo) + -q] . (4) 

The real valued parameters determining this 'gap soliton' 
are $7 (frequency), (j) (initial phase) and zq (center), and 
we have 



= cos q 



(5) 



Then A plays the role of the 'wave number' oi the evanes- 
cent wave. 

Nonlinear supratransmission threshold. For a Bragg 
medium, extending in the region z e [0,L], and initially 
in the dark, we set the initial data 



e(z,0) = 0, /(z,0) = 



(6) 



The boundary value problem that mimics the scattering 
of an incident radiation at frequency loq + kcQ on the 
medium in z = is 



e{0,t) ^ Ae 



-int 



fiL,t)^0, 



(7) 



where the second requirement means no backward wave 
incident from the right in z = L. The constant A (in 
general complex valued) is the dimensionless amplitude 
of the incoming radiation. Equations © and {Tj) consti- 
tute a well posed initial-boundary value problem for the 
PDE @. 

In the linear case such a boundary forcing would simply 
generate the evanescent wave e{z, t) = A g^*^*^-*^^ . In the 
nonlinear case however this boundary forcing generates 
the solution Q for the value of (f> and zq < such that 
the amplitude in z = be precisely A. Then, for each 
fixed forcing frequency f2, there exists a maximum value 
As of A beyond which there is no solution {0, zq}. This 
threshold As is given by |e(0,i)| from Q evaluated at 
Zq = 0, namely by 



A, = 2x/273 sin(- 



'2/3 sin(- arccosJl) 



(8) 



This is the threshold amplitude of the incident envelope 
above which the system develops an instability and gen- 
erates a propagating nonlinear mode. 




FIG. 1: Plots of the amplitude |e(2,t)| of the right-going en- 
velope at given time (t = 145) for two values of the driving 
amplitude A at frequency f2 = 0.8. 



Nonlinear dispersion relation prediction. By seeking 
a solution of (01 e = aexp[i{kz — cot)], f = 6exp[i(fcz — 
ujt)] one gets a nonlinear algebraic system for the un- 
knowns a (incident amplitude) and lo (incident fre- 
quency) expressed in terms of the wave number k once 
the amplitude b has been eliminated. This gives the non- 
linear dispersion relation uj{a, k) as a solution of a third 
order algebraic equation. One of the 3 solutions must be 
discarded as it is singular in the linear limit a — s- 0, the 
other two giving the deformations of the linear branches 
Lu{k) = ±Vl + k^. It can be shown that they correspond 
to the relations b = —a (upper branch) and b ~ a (lower 
branch) . 

Being interested in the shift of the gap, it is sufficient to 
study the solutions at fc = where the system simplifies, 
and to stick with the upper branch (as we look at right- 
going incident waves) for which b — —a. Then ^ readily 
provides the value w(0, a) = 1 — | of the gap opening. 

For our purpose it is more convenient to express the 
value Ara of the incident amplitude a for which the inci- 
dent frequency touches the nonlinear gap edge w(0, a), 
namely 



A„, = \j-{l-n) 



(9) 



This formula provides a prediction of transparency by 
nonlinear shift of the gap. It is now compared to our 
prediction (jS)) by means of numerical simulations. 

Numerical simulations. The system ^ is solved by 
means of an semi-implicit third order finite difference 
scheme and boundary values are taken into account at 
first order. Explicitely we first rewrite Q for two func- 
tions u(z, t) and w(z, t) and replace the operators dt + dz 
by the set of differences for M„(t) = u{nh,t) 

1/ ^ • 1 . 

Ui + -(Ui~ Uq), U2 + —{U2- Uq), ■■■ 

h 2n 



1 



1 / ^ • 1. 

""AT-l + 7rr("JV+l - Un~i), UJV + T{''^N " UAT-l), 

2n h 

so as for f„(i) = v{nh,t). The length is L ~ hN and 
overdot means time differentiation. Equation |2Jl results 
as a system of 2N coupled ODE then solved through the 
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subroutine dsolve of the MAPLES software package that 
uses a Fehlberg fourth/fifth order Runge-Kutta method. 
Finally the solution, e.g. e{z, t), is obtained from u„{t) by 
e{z,t) = [un+i{t) + Un{t)]/2. Although such code sends 
some numerical noise in the solution, it is uneffective for 
sufficiently small h depending on the required time of 
integration. 

Next, in order to avoid initial shock, the boundary con- 
dition e(0, t) is smoothly turned on and smoothly turned 
off by assuming instead of ||7)| 



mitted (T) energies at given arbitrary time tm by 



e{0,t) = — [tanh(p(i 



to)) - tanh{p{t - ti))] . (10) 



A practical interest of such an incident wave is that it 
reproduces the nuasi- constant wave irradiation of the ex- 
periments of |13| . 

Most of the results presented here are obtained with 
N = 120 spatial mesh points, h = 0.05 grid spacing 
over (hence a length L = 6 normalize units) a time of 
integration tm = 200. The parameters of the boundary 
field are a real-valued amplitude A, a slope p = 0.2, an 
ignition time to — 20 and an extiction time ti = 180. 

We display in figure ^ an instance of two different nu- 
merical solutions (the modulus of the right-going enve- 
lope) for an incident frequency = 0.8 and amplitudes 
A = 0.51 (no transmission) and A = 0.52 (gap soliton 
generation) when the theoretical threshold predicted by 
© is ^ = 0.5164. 

Using this simple diagnostic, the bifurcation predicted 
by expression ijSJ) is numerically checked for a series of 
frequency values (in the range [0.1, 0.995]) and we obtain 
the figure |2] where the dots represent the smallest value 
(with absolute precision of 10^^) of the amplitude A for 
which nonlinear supratransmission is seen to occur. The 
expression © is also plotted (dashed line) which shows 
that the nonlinear shift provides a wrong answer. 




FIG. 2: Values of the threshold of nonlinear supratransmis- 
sion observed on numerical simulations (crosses) as compared 
with expression (|HJ (solid line) and with @ (dashed line). 

Bifurcation of transmitted energy. The system ^ 
possess the conservation law 



5t(|eP + |/n + 9.(|ep-|/P) = 0. 



(11) 



For a given boundary condition ||7J) (i.e. for fixed A and 
O), we define the incident (/), reflected (R) and trans- 



i{A,n) = 
R{A,n) = 



dt\e{0,t)\^ 



dt|/(0,t)p 



T{A,n) = / dt |e(L,t)p 



(12) 

(13) 
(14) 



Then, upon integration of Ijllf) on the length [0, L] and 
time [0, tfn] we obtain 



R + T-I+ / dz (|e(z,t™)|2 + |/(z,i„)|2) = . (15) 
Jo 

If tm is larger than the irriadiation duration, the en- 
ergy injected eventually radiates out completely, namely 
e{z,tm) and f{z,tm) vanish, and we are left with the 
conservation relation R + T = I which can be written 



p(A,r!)+r(A,r!) = 1, p-i?//. 



T/I. (16) 



The figure 13 shows a typical numerical simulation 
where, for = 0.95, we have computed the reflection and 
transmission factors, together with their sum, for 50 dif- 
ferent values of the amplitude A in the range [0.20, 0.45]. 
These simulations showt the sudden energy flow through 
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FIG. 3: Plot of the reflection factor p (full line), transmission 
r (dashed line) and sum p + r (dotted line) for system Q 
with boundary value IIUII and parameter values given there. 

the medium, as a result of nonlinear supratransmission. 
Note from the sum p + r before the bifurcation on figure 
13 that the numerical code conserves the energy with a 
precision of 2%. 

The transmissivity is due to the generation and prop- 
agation of light pulses shown figure ^ representing the 
energy density ]e(L,t)]^ measured at the output as a 
function of time. We have also plotted the input energy 
density given in H10|l with parameters t^ = 100, to = 20, 
ti = 90 and p = 2. We show now that these are gap 
solitons travelling at a fraction of the light velocity (slow 
light pulses). 

Travelling gap solitons. Although the system (1^ is 
not Lorentz invariant, the propagating solution of 
can still be written in terms of the boosted variables 119| ■ 
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FIG. 4: Energy density output \e{L,t)\^ (full line) compared 
to the input (dashed line) for lj = 0.95 and A = 0.5. 



As we are intersted here only in the energy flux, we write 
the soliton solution |e(2;,i:)p moving at velocity v < 1 
with frequency cosq as 



\e{z,t)\' 



1/2 



cosh 



(17) 



Such an expression allows to fit a given simulation by 
seeking the two parameters v and q, and the initial po- 
sition zq, that reproduce with the explicit solution the 
results of the numerical simulation. 




FIG. 5: Fit of the numerical solution (crosses) by analytical 
soliton solution (solid line) at times 65.7 and 72.7, for an 
incident frequency fl = 0.8 and amplitude A — 0.53. 



An instance of such a fit is displayed in figure O that 
shows the analytic soliton |e(2:,t)p compared with the 
numerical simulation for two fixed values of time. This 
furnishes the following velocities 
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These are slow light pulses and naturally, small velocities 
are obtained close to the gap edge where the required 
energy is small. 

Conclusion. The property of a Bragg medium in the 
nonlinear Kerr regime to sustain nonlinear supratrans- 
mission has allowed us to obtain the analytic expression 
^ that fixes the threshold amplitude of an incident cw- 
beam in terms of its departure from the Bragg resonance. 

This constitutes a practical tool to investigate switch- 
ing properties of a Bragg medium and allows at the same 
time to understand the process at the origin of sudden 
transmissivity of the Bragg mirror. It is the nonlinear 
supratansmission that results from a nonlinear instabil- 
ity intrinsic to boundary value problems jl^ . 

We expect our result to be useful! for experiments by 
the ability of the Bragg medium to become partly trans- 
parent by gap soliton generation for an incident beam of 
low energy if its frequency is chosen close to the gap edge. 
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